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We give a proof of the KAM theorem on the existence of in- 
variant tori for weakly perturbed Hamiltonian systems, based 
on Thirring's approach for Hamiltonians that are quadratic in 
the action variables. The main point of this approach is that 
the iteration of canonical transformations on which the proof 
is based stays within the space of quadratic Hamiltonians. 
We show that Thirring's proof for nondegenerate Hamiltoni- 
ans can be adapted to Hamiltonians with degenerate twist. 
This case, in fact, drastically simplifies Thirring's proof. 

PACS numbers: 03.20.+i, 05.45. +b 



I. INTRODUCTION 

Regular Hamiltonian dynamics is characterized by the 
existence of as many independent conserved quantities as 
degrees of freedom d. As a consequence, each trajectory 
is confined to evolve on an invariant torus of dimension d. 
The KAM technique was developed to prove the stability 
under small perturbations of a large fraction of these in- 
variant tori. The KAM theorem states that, if the 
frequency satisfies a diophantinc condition, and the size 
of the perturbation e is sufficiently small, then a torus of 
that frequency will be stable. The proof is based on an 
iterative algorithm to construct the invariant tori. Each 
step of the KAM iteration consists of a coordinate trans- 
formation that reduces the size of the perturbation from 
order e to e 2 . 

In this paper, we derive a KAM theorem for a family of 
Hamiltonians that has been used in Refs. within 
the setup of a renormalization-group approach to the 
breakup of invariant tori. 

We consider the following class of Hamiltonians with d 
degrees of freedom, that are quadratic in the action vari- 
ables A = (Ai, A2, . • • , Ad) £ R d and described by three 
scalar functions of the angles <p — (<pi,<p2, ■ ■ ■ , <Pd) £ T d 
(the d-dimensional torus parametrized, e. g. by [0, 2n] d ): 
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No 



^r^K, V I / = ("i,.,^)e2 i \{0}, (1.2) 



ff(A, 9 ) = ^)(!1.4) 2 



[uj +g(<p)fl}- A + f(<p) 



(1.1) 



where ujq 6 R d is the frequency vector of the considered 
torus, and fl is some other constant vector. Without 
loss of generality, we assume that fl = (fix, . . . , fid) is of 

norm one, i. e. \ fl\ = J2i=i l^il = 1- The functions m, g, 
and / are real analytic on T d , i. e. they are holomorphic 
functions on some complex neighborhood of T d . 
We consider Wq verifying a diophantinc condition 



for some r > d — 1 and a > 0. 

The aim is to prove the existence of a torus with fre- 
quency vector <jJq for Hamiltonian systems described by 
Eq. (1.1). The method is to find a canonical transfor- 



mation such that the equations of motion expressed in 
the new coordinates show trivial ly t he existence of this 
torus. For Hamiltonians of type (1.1), if one takes g and 
/ equal to zero, then the resulting equations of motion 
are 



dA _ 1dm 2 
dt ~ 2 8<p [ ' ' 
dip 

-Zj- = m(ip)(fl ■ A)fl + uj . 



(1.3) 
(1.4) 



Thus, there exists a torus with frequency vector u> lo- 
cated at A — (even if the resulting Hamiltonian is not 
globally integrable, i. e. for A ^ Q). This canonical 
transformation cannot be defined directly, because the 
formal expressions that appear in the classical perturba- 
tion theory do not converge due to the presence of small 
denominators. The construction is done via an iterative 
algorithm. We iterate a canonical chan ge o f coordinates 
that maps a Hamiltonian of the form (1.1) with a per- 
turbation (g,f) of order 0(e), into a Hamiltonian with 
(g 1 , /') of order 0(e 2 ). When the iteration converges, the 
perturbation (g, /) is completely eliminated. 
The fact that m does not need to be eliminated to prove 
the existence of the torus with frequency vector u)q, al- 
lows us to stay with Hamiltonians that are quadratic 
in the actions at each step of the iteration. This ap- 
proach has been developed by Thirring In fact, 
Thirring considered a non-degenerate family of Hamilto- 
nians quadratic in the actions, of the form 



H(A,<p) = -A-M(<p)A 



> + g(tp)] -A + f(<p) 



(1.5) 



where M is a d x d matrix such that det M / 0, a nd g 
a vector. This implies that Thirring's Hamiltonian (LL5J ) 
satisfies the "twist condition" 



det 



d 2 H 



dAdA 



det M ^ 0. 



(1.6) 



The Hamiltonian ([bl]) does not satisfy the twist con- 
dition for any (p. The extension of the KAM theorem 
for de gen erate systems was done by Arnold |l(| (see also 
Refs. ]TT|-^5[). In the case we consider here, the rank 



1 



of the matrix in (1.6) is one. Thus, there is a twist in 
a particular direction, which is the only relevant one for 
the conside red perturbation. 

Condition (1.6) is also called (standard) nondegeneracy 
condition. KAM theorems were also proven for Hamilto- 
nians which satisfy the isoenergetic nondegeneracy condi- 
tion jl6|,[l7| , for which the following determinant of order 
d + 1 docs not vanish 



det 



d 2 H OH 
dAdA dA 

OH 
~dA 



(1.7) 



For d — 2, the determinant ( |l.7| ) is equal to 
— m(<p)[det (n,uj )] 2 . Thus, the isoenergetic nondegen- 
eracy condition is not sat isfied if m(ip) has zeroes. 
For d > 2, the model (IT) is isoenergetically degenerate. 
The present result is thus not a direct consequence of the 
ones in . We notice that in t he case where fl is par- 

allel to u>q, the Hamiltonian (IT) is integrable (because 
of the existence of d integrals of motion in involution : 
H and u}q ■ ip, where u3q denotes a vector perpendicular 
to u) , and there are d — 1 such independant vectors). 
In this article, we present a sel f-co ntained proof of the 
KAM theorem for Hamiltonian (1.1) based on Thirring's 
proof for Hamiltonian (fh^). The advantages are twofold: 
On one hand, we show that Thirring's proof can be 
adapted to degenerate twist Hamiltonians, and on the 
other hand, the resulting proof becomes even simpler 
than Thirring's. The fact that the iteration stays within 
the space of Hamiltonians quadratic in the actions is very 
useful, e. g., for numerical implementations to study the 
breakup of invariant tori |5|-|7]] . 

Before entering into details of the theorem, we give ba- 
sic definitions and notations: We denote T) p a complex 
neighborhood of T d defined by 



V p = { V eC d \ ||Im^||<p}, 



(1.8) 



where ||z|| = maxi(|zi|), for any z € C d . We will consider, 
in the following calculations, scalar functions defined in 
T> p . More precisely, we define A p as the set of complex 
functions f(<p) defined on T> p , analytic in the interior 
of T> p , of period 2n in the variables tpi, and which have 
real values when ip e R d . We define a norm on A p : 
II/IIp = sup^g-p \f(<p)\. We define (/}, the mean value of 
/by 



d d p 
(2^p 



/(¥>)■ 



(1.9) 



In the following sections, we will use the notation df = 
df 

— — for any function of the angles. 
dp 

In Sec. y, we define the KAM iteration. In Sec. |HJ, 
we give estimates on the transformed functions. Finally, 
in Sec. GVL we iterate the transformation, and prove the 



f ollo wing KAM theorem for the family of Hamiltonians 



Theorem 1 For H(A,<p) = -m(<p)(Sl ■ A) 2 + [w + 

g(p)Q] ■ A + f(<p), suppose that 

\i) |u; • u^ 1 < a\u\ T , Vi/ G Z d \ {0}, for some a > 
and t > d — 1; 

(ii) m, g, and f are analytic in D p , and |0| = 1; 

(Hi) max(||g|| p , ||/|| p ) < f3 , where 

A) = 2- 10 r- 6 c- 3 (V3) 3(r+d+1) , 

c = 2 3d a [2e^(r + l)] T+1 , 
T = max(l, ||m|| p , Km)! -1 ), 



h < inf(2p/9, 



J/(r+d+l) V 



Then, there exists a canonical transformation, analytic 
in T> p _<j h / 2 , such that the Hamiltonian expressed in the 
new coordinates is 

= ^m(°°\<p^)(n ■ A (oo) ) 2 + w • A iao \ 

where m'°°' is analytic in T> p _ 9h /2- As a consequence, 
the system has an analytic invariant torus of frequency 

LJ . 

Remark: if the constant part (m) of the quadratic term 
is zero, the maximal amplitude of the perturbation [3q 
given by the theorem becomes zero. Thus, in order to 
have a nontrivial result, we require that (m) is nonzero. 



II. EXPRESSIONS OF THE GENERATING 
FUNCTION AND OF THE NEW HAMILTONIAN 

We perform a canonical transformation Uf ■ Up, A) i— > 
(ip',A') defined by a generating function [|l8|,[l9| linear 
in the action variables, and characterized by two scalar 
functions Y, Z, of the angles, and a constant a, of the 
form 

F(A', <p) = (A' + aff) ■ <p + Y(tp)fl ■ A' + Z(<p), (2.1) 
leading to 

dF 

A= — =A' + (n-A')dY + a n + dZ, (2.2) 
o<p 



dF 
d~A i 



ip + Y(tp)n. 



(2.3) 



Inserting Eq. ( |2.2| ) into the Hamiltonian ([O]), we obtain 
the expression of the Hamiltonian in the mixed represen- 
tation of new action variables and old angle variables 



H{A' , tp) = -m(tp)(tt ■ A') 2 



-[cj + g{<p)Sl]-A' + f{(p) 



(2.4) 



2 



with 



rh = + dY) 2 m, 

g = g + uj Q ■ dY + mb + £1 ■ dY (g + mb) , 



f = f + ujQ-dZ- 



1 



i& 2 



gb, 



(2.5) 
(2.6) 

(2.7) 



where b(<p) = att 2 + fl ■ dZ. As the new an gles do not 
depend on the actions, the Hamiltonian (1.1) expressed 
in the new variables is also quadratic in the actions, and 
of the same form as (1.1). We notice that this transfor- 



mation does not change fl. 

We determine the generating function (2.1) such that the 



new functions g, /inHo Uf vanish to the first order in 
E. This leads to the conditions 



u)q ■ dZ + / = const, 

ujQ-dY + g + m (an 2 + n ■ dZ) 



0. 



(2.8) 
(2.9) 



The constant a allows us to have (<?} = 0(e 2 ), in order 
to keep the frequency u> at the chosen value. We recall 
that the functions g and / are of order 0(e) and m is of 
order one; as a conse quen ce Y, Z and a are of order 0(e). 
Equations (2.8) and (2.9) are solved by representing them 
in Fourier space. They define the generating function F 
as 



uj ■ v 



(g) + (mfi • dZ) 
Y ( ( P) = £ ^—(9v + (mfi • 8Z) V 



(2.10) 



(2.11) 



+TO„af2 2 ) 



2\ e iv-ip 



(2.12) 



where the scalar functions m, <?, f are represented by 
their Fourier series, e. g. 



Thus the scalar functions in H become 
m = (1 + n ■ dYfm, 



g = - (w • BY) (n • dY) , 

f = \{g-uo-8Y) (an 2 + n-dz) 



(2.13) 



(2.14) 
(2.15) 

(2.16) 



The expression of H in the new angles requires the in- 
version of Eq. (2.3). We denote H' the Hamiltonian ex- 
pressed in the new variables H'(A' , tp') = H(A' , ip): 

H>(A'^')= l - m \ v '){n-A') 2 

+[u +g'(<p')n]-A' + f(cp') . (2.17) 

In the following section, we give estimates on Z , Y and 
a, in order to derive estimates on m', g' and /'. 



III. ESTIMATES ON THE GENERATING 
FUNCTION AND ON THE NEW HAMILTONIAN 

A. Estimate on the generating function 

Lem ma 1 Let f be element o f A B , and ujq satisfy Eq. 

then Z((p) given by Eq. (2AL ) is analytic on T> p -h 
(for any choice of h with < h < p < 2), and satisfies 
the following estimates 



wht 



\\z\\ p - h < ch-^WfH, 

\\u> -dZ\\ p - h < \\f\\ p , 

\\n-Bz\\ p - h <ch- T - d - 1 \\f\\ p , 

2 3d a (2e- 1 r) T and c = 2 3d a ^e" 1 ^ + 1)] ' 



(3.1) 
(3.2) 
(3.3) 



Proof: 

This lemma is proved, for instance, in Refs. ]To| , pfH . 
First, we estimate the Fourier coefficients of / expressed 
as integrals of / over the torus, 



d d ip 



(3.4) 



By a shift of the angles ip t— > tp+irj, where r\i = —pvi/\i/i\, 
d d ip 



fu 



d (2tt) 



-/( ¥ , + i 77 )e— <P e 



iv-tp —p\u\ 



(3.5) 



Then |/,| < ||/|| p e-^l-l, for all ^ e 

To estimate Z given by Eq. (2.10), we use the diophantine 
property (|J), 



\Z(<p)\<aJ2W\ T \Me- u - Imv . 

Then, for all <p £ 2?p-/i, we have 

\Z( V )\<a\\f\\ p J2H T e- hH - 



(3.6) 



(3.7) 



To estimate the sum, we use the following property which 
is easy to check (see Ref. 

\v\ T < (^) e^ h ' 2 , W,h>0. (3.8) 

From the fact that 1/(1 - er x ) < 2/x, for all x G]0, 1[, 
we have 



(3.9) 



which gives the estimate (^lj). The estimate ( |3.3j ) is ob- 
tained by the same calculati ons , and the estimate (3.2) 
is straightforward from Eq. (2.8). □ 
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Lemma 2 Let m, g, and f be elements o f A , and ujq 
satisfy Eq. (Li), then Y(ip) given by Eq. ( 2.1i ) is ana- 
lytic on T> p _2h- The constant a and the function Y satisfy 
the following estimates 

\\Y\\ p . 2h < ch-r- d (l + Km>r x ||m|| p ) 

x(||$L + c/r T - d -i/|| p |M| p ), (3.10) 

||«0 • 0^11,-2* <(l+|(m>rVllp) 

xdlflllp + cft-^-^H/llplHIp), (3.11) 



the width of the strip T> from p — 2h to p — 3ft, e. g. 
||/'[[/)-3/i < ||/||p-2/i- Then, the estimates on t he new 
perturbation (<?', /') are obtained from Eqs. ( 2.15 )-( 2.1€ ) 
and from the estimates on the generating function: 



\\g'\\ P -sh < 2Wh-^ +d+1 \ 2 , 



(3.19) 
(3.20) 



\Sl-dY\\ p - 2h < ch 



-T-d-l 



l + IWi^iHip) 



x (Il5llp + cft 



-T-d-li 



Hip), (3.12) 



< Km))- 1 (\\g\\ p + eft— ||m|| p ) . (3.13) 



The proof is the same as the one for the estimates on Z. 



B. Estimate on the new Hamiltonian 

The expression of the scalar functions of the Hamilto- 
nian expressed in the new actions and old angles are ex- 
plicitly known, and estimates on these functions are easy 
to obtain from the estimates on the generating function. 
One has then to invert Eq. (2.3) in order to obtain the 
estimate on the Hamitonian expressed in the new angle 
variables. The Jacobian of this transformation is 



dct 



1 + il-dY. 



In this section, we denote V p — max 
max(l, ||m || , Km)! -1 



(3.14) 



\g\\ p ) and V 



From Eq. (2.3) and estimate (3.10), one has the following 
inequality: 



\cp' - (p\ < 2Y 3 ch 



-r—d , 



1 + ch 



-r-d-l\ 



(3.15) 



for all <~p 6 T> p -2h (we recall that Y > 1). If we assume 
that V p and h are sufficiently small, one has an estimate 
on the new angles. More precisely, we assume the follow- 
ing inequalities: 



T 3 c 2 h -2(r+ d +l) v < 



1 



eft 



-T-d-l 



> 1. 



Then, for all ip G 2? p -2h, we have 



(3.16) 
(3.17) 

(3.18) 



As a consequence, T> p _ 2 h C $(D p _ 3 h) by the map <p' i— > 
tp = &(>p>') given by Eq. fl2.3| ). In order to express the 
estimates with respect to ^7 it suffices thus to restrict 



Taking into account condition (3.17), we obtain the esti 
mate on V' = : 



p—Zh 



(llff'llp-aOTVsh): 



V' p _ 3h < 2 4 Y 6 c 3 h-^ T+d+1 W p 2 . 



(3.21) 



Concerning the quadratic term, we deduce from Eqs. 
(§T2|) and that 



\\m'\\ p - 3h < \\m\\p{l + AY 3 c 2 h- 2 ^+ d +^Vpy . (3.22) 
Then using Hypothesis (Hi) of the theorem, we obtain 

IKII„-3h < 2r. 

The mean value of m! is determined by the integral 



<m') = 



d d ip' 
Wi' 



yd 



■m'(<p>). 



(3.23) 



With the change of v ariabl e <p' t— >• tp given by Eqs. (|2 
( 3.14 ) and using Eq. ( 2.14 ), we rewrite the integral 



(m') = 



d d cp 
,. d (2ir) d 
d d ip 



(2n) 



1 + fl ■ dY\fh(ip) 



+ dY\ 3 m(tp). 



(3.24) 



To estimate \(m')\ 1 , we fi rst estimate the difference 
l(m') — (m)| using Eq. ( 3.24 ), and we obtain 



l(m') - <m)| < ||m||p(3||n • dY\\^ 2h + 3||0 • dY\\ 2 p _ 2h 
+ \\n-dYf p _ 2h ). (3.25) 



From Eq. (3.12) and condition ( 3.16| ), we have the esti- 
mate 

||n • dY\\ p _ 2h < 4r 3 c 2 ft- 2(T+<i+1 V p < 1, (3.26) 

which leads to 

|(m') - (m)| < 28r 4 c 2 ^ 2(T+d+1 V p . (3.27) 

One can easily check that from Hypothesis (Hi) , it follows 
that 



|(m') - (m)| < 



1 

2T ' 



(3.28) 



Writing that (to') = (to) + (m') — (to), we have the fol- 
lowing estimate on |(m')| _1 : 



IK)!" 1 < 



Km)!- 1 



(3.29) 



|l-|<m)|-i|(m')-(m)|f 

From Eq. ( |3.28[ ), we deduce that Km')!" 1 < 
2r. Therefore, we have proved that Y' = 
max(l,||m'||p_ 3 ^,Km')|- 1 ) < 2Y. 
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IV. CONVERGENCE OF THE ITERATION: KAM 
THEOREM 



The estimate ( 3.21 ) gives a precise meaning to the 
statement that the perturbation is reduced from V to 
V 2 : we see that the actual reduction is somewh at sm aller, 
since h depends on V p through the condition ( 3.1 6j ) [the 
actual rate of reduction that takes this into account will 



be expressed by Eq. (4.11) below]. The counterpart is 
that the width of the domain of analyticity of the new 
functions is reduced from p to p — 3h; for this reason, h is 
called the "analyticity loss parameter" . At the nth step 
of the iteration, we choose h n as the analyticity loss pa- 
rameter such that the final domain of analyticity (after 
an infinite number of iterations) does not shrink to zero. 
For instance, we choose 



hn = h3> 



with h < inf (2/9/9, c 



l/(r+d+l)N 



(4.1) 



The second term in the inf is to guarantee also the con- 
dition 

( p7| ). We denote T n = max(l, ||mW|| Pn , Km^)! -1 ) 
and V n = max(||£/W|| Pn) ||/ ( " ) || P „), where mS n \ g {n) and 
de note the three scalar functions defining the Hamil- 
tonian (1.1) after n iterations, and 



Pn-i — 3h n 



3£> 

k=l 



(4.2) 



The previous section gave the following estimates 

V n <2^ n _^h- 3 ^^_ 1 , (4.3) 

||m(")|U < Hm^- 1 ) |U_ t 

x(l + 4rf l _ lC 2 / l - 2 (-+ d + 1 )K-i) 2 . (4.4) 

We now prove that T n is bounded, and that V n tends to 
zero faster than geometrically as n goes to infinity. 
Denoting 



V = maxdlfft, \\f\\ p ), 

T = max(l, \\m\\ p , |(m)| _1 ), 

7 = 2 10 r 6 c 3 /i~ 3(r+d+1) , 

6 = 3 3(x+d+l) ) 

we define the sequence of values 

(jSVo) 2 ' 1 



We will show by induction that 

r fc < 2r, 

Vfc < £fe, 



7<5 



n+l 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



(4.9) 



(4.10) 
(4.11) 



for all k. If we assume Vb to be sufficiently small, such 
that 76V0 < 1, i. e. 



2 10 r 6 c 3^-3(r+ ( i+l)33(r+d+l) yo < ^ 



(4.12) 



which is the hypothesis (Hi) of the theorem, then V n 
tends to zero faster th an S ~ n = 3 ~ 3 "( T + rf + 1 ) as n goes to 
infinity. The bounds (fh^- ([□!]) are satisfied for k = 0, 
since T = T > and Vb = £o- 

Suppose th at f or all k < re, these bounds are satisfied. 
From Eqs. (O) and ( 4.10 ), we deduce that 



(4.13) 



The n usi ng Eq. ( 4.11 ), we have V n < e n - Concerning T„, 
Eq. (O) leads to 



l|m (n) IL < l|m| 



p 

n-l 



Y[ (l + 2 5 T 3 c 2 h~ 2{T+d+1) 3 2k(T+d+1) V k 



k=0 



Using Eq. (4.11) and Hypothesis (in), we obtain 
2 5 T 3 c 2 hr 2< - T+d+1 ' ) 3 2k ^ T+d+1 ' > Vk < 3~( k + 3 )( T + d + 1 ) . (4.14) 
Then we have 

2 



\m 



(«)| 



< 



r y[ (1 + 3-( fc+3 H T + d+i )y 



(4.15) 



k=0 



Using the fact that J|(l + Xk) 2 < exp2J2 x k, we show 
that the infinite product converges, and that it is smaller 
than 2. Thus, Hm^H^ < 2r. 
For |(m' n ))| — x , we use estimate (3.29): 



< 



1 



Km*"*)!" 1 

|( m (n-l))|-l - |1 - |( m («-l))|-l|( m (n)) _ ( m (n-l))|| 

(4.16) 

The differ ence \{m^ n ') — (m(" _1 ))| is evaluated as in the 
first step ( ^27| ): 



(n) 



) - (mC™- 1 ))! < 28 • 2 4 T 4 c 2 h- 2 ^ +d+1 ^ 

x3 2(n-l)(r+«tfl)^ n _ 1> (417) 



Using Eq. ( 4.11 ), Hypothesis (Hi), and the fact that T 
and ch~( r+d+1 ) are greater than one, 

|( m W) - (m^- 1 ^ < 2- 1 r- 1 3- (n+2)(T+d+1) . (4.18) 

Then, 



|(m("))r 1 <rn(l-3- 



(fc+2)(T+d+l) 



k=0 



(4.19) 



Using the fact that J|(l — x^) -1 = exp(— ^ln(l — x/.)), 
and that - ln(l - x) < x2 In 2, Vx £ [0,1 /2] , we show by 
straightforward calcul ation s that | (to'™)) | _1 < 2r. 
Thus we have proved (4.10)-(4.11) for all k. 



We finally have to verify that with the choices ( 4.1 )-( 4.2 ) 
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the conditions (3.16)-( p.l7[) a re satisfied at each st ep o f 
the iteration. Equation fl3.17|) is guaranteed by Eq. (4.1): 

chn {T+d+1) > c/i-( T+d+1 ) > 1. For Eq. (fl(}), we first 



notice that, since T > 1, T n < 2T, and chr^ r+d+1 ^ > 1, 
we can write 



Therefore, using Eq. (4.11) and Hypothesis (Hi), we have 

1 



< 



We have shown that, for sufficiently small initial pertur- 
bation, the Hamiltonian converges under successive iter- 
ations of the transformation to a Hamiltonian of the form 
(1.1) with g = f = 0. In order to complete the proof, 
we have still to verify that the composition of the in- 
finitely many canonical transformations is a well-defined 
finite canonical transformation. It suffices to verify that 
\ip(°°) - p\ an d \a(°°* > - A\ are finite. This is an immedi- 
ate consequence of the fast convergence of the iteration. 
From Eq. (3.15), we deduce e. g. that 



<^2h n = 3h/2. 



(4.20) 



n=l 



Analogously, using Eqs. ( p^ ), ( ft.lip , ( |3.12j ), and 

(4.11), for \A\ < A , we can bound 



|A (oo) -A\< const x ]T h n . (4.21) 
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